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Abstract. Grapii transformation is concerned with the manipulation of graphs 
by means of rules. Graph grammars have been traditionally studied using tech- 
niques from category theory. In previous works, we introduced Matrix Graph 
Grammars (MGGs) as a purely algebraic approach for the study of graph gram- 
mars and graph dynamics, based on the representation of graphs by means of 
their adjacency matrices. MGGs have been succesfuUy applied to problems such 
as applicability of rule sequences, sequentialization and reachability, providing 
new analysis techniques and generalizing and improving previous results. 
Our next objective is to generalize MGGs in order to approach computational 
complexity theory and static properties of graphs out of the dynamics of certain 
grammars. In the present work, we start building bridges between MGGs and 
complexity by introducing what we call Monotone Complex Logic, which allows 
establishing a (bijective) link between MGGs and complex analysis. We use this 
logic to recast the formulation and basic building blocks of MGGs as more proper 
geometric and analytic concepts (scalar products, norms, distances). MGG rules 
can also be interpreted - via operators - as complex numbers. Interestingly, the 
subset they define can be characterized as the Sierpinski gasket. 



1 Introduction 

Graph transformation ||T6| is concerned with the manipulation of graphs by means of 
rules. Similar to Chomsky grammars for strings, a graph grammar is made of a set of 
rules, each having a left and a right hand side (LHS and RHS) graphs and an initial host 
graph, to which rules are applied. The application of a rule to a host graph is called a 
derivation step and involves the deletion and addition of nodes and edges according to 
the rule specification. Roughly, when an occurrence of the rule's LHS is found in the 
graph, then it can be replaced by the RHS. Graph transformation has been successfully 
applied in many areas of computer science, for example, to express the valid structure 
of graphical languages, for the specification of system behaviour, visual programming, 
visual simulation, picture processing and model transformation (see (31). In particular, 
graph grammars have been used to specify computations on graphs, as well as to de- 
fine graph languages (i.e. sets of graphs with certain properties), thus being possible 
to "translate" static properties of graphs such as coloring into equivalent properties of 
dynamical systems (grammars). 

In previous work II12I13I14I15I we developed a new approach to the transformation 
of simple digraphs. Simple graphs and rules can be represented with Boolean matrices 



and vectors and the rewriting can be expressed using Boolean operators only. One im- 
portant point of MGGs is that, as a difference from other approaches |fT6l . it explicitly 
represents the rule dynamics (addition and deletion of elements), instead of only the 
static parts (pre- and post- conditions). Apart from the practical implications, this fact 
facilitates new theoretical analysis techniques such as for example checking indepen- 
dence of a sequence of arbitrary length and a permutation of it, or obtaining the smallest 
graph able to fire a sequence. See I.15J for a detailed account. 

In lfT4l we improved our framework with the introduction of the nihilation matrix, 
which makes explicit some implicit information in rules: elements that, if present in 
the host graph, disable a transformation step. These are all edges not included in the 
left hand side, adjacent to nodes deleted by the rule (which would become dangling) 
and edges that are added by the production, as in simple digraphs parallel edges are 
forbidden. In this paper, we further develop this idea, as it is natural to consider that a 
production transforms pairs of graphs, a "positive" one with elements that must exist 
(identified by the LHS), and a "negative" one, with forbidden elements (identified by 
the nihilation matrix). 

Complexity theory [611 H is concerned with the study of the intrinsic complexity of 
computational tasks. Traditionally, it has been studied through abstract devices able to 
represent the notion of algorithm, such as Turing Machines or Boolean Circuits f\n\. 
Our proposal is to use MGGs instead, as its algebraic nature allows using results from 
different branches of mathematics such as logics, group theory and Boolean algebra. 

In this paper we give a first step in the direction of approaching complexity the- 
ory with MGGs, by introducing Monotone Complex Logic (MCL). Similar to complex 
numbers, a complex formula in MCL contains a certainty and a nihil part, both Boolean 
propositional formulas. We use MCL terms to encode the "positive" part of a simple 
digraph (the LHS) and the elements that cannot be found (e.g. the nihilation matrix). 
Using a rational encoding of adjacency matrices, we can express complex terms re- 
ferring to simple digraphs into the unit interval of complex numbers C([0, 1]). Interest- 
ingly, the set of complex numbers defined by valid MCL terms on simple digraphs is the 
well-known Sierpinski gasket fractal JS). The rational encoding allows using geometric 
and analytic concepts, for example, we have defined a xor-based norm for MCL terms 
which can be interpreted as the number of elementary operations needed to transform 
one digraph into another 

Thus, we can use MCL terms to redefine and extend all concepts of MGGs. In 
this paper we introduce the encoding of productions in both its static and dynamic 
formulations. In the dynamic formulation of a production, the rule dynamics (element 
addition and deletion) are also represented as an MCL term, and thus belong to the 
Sierpinski gasket too. We also show the generalization of the main MGG concepts, like 
coherence, compatibility, initial digraphs, image of sequences and G-congruence using 
MCL. 

Paper organization. Section |2] gives a brief overview of the basic concepts of MGGs. 
Section|3]introduces MCL, used to establish a link between MGGs and complex analy- 
sis. Section|4]encodes graphs as complex numbers, and a scalar product, a norm and a 
notion of distance are introduced. Section |5] encodes productions as complex numbers 
and completes the link between MGGs and complex numbers. Sections |6] and [T] gen- 



eralize the main sequential results of MGGs such as coherence, compatibility, initial 
digraphs, image of sequences and G-congruence. Finally, Sec.|8]ends with the conclu- 
sions and further research. 

2 Matrix Graph Grammars: Basic Concepts 

In this section we give a very brief overview of some of the basics of MGGs, for a 
detailed account and accesible presentation, the reader is referred to [15 J. 
Graphs and Rules. We work with simple digraphs, which we represent as (M, V) 
where M is a Boolean matrix for edges (the graph adjacency matrix) and V a Boolean 
vector for vertices or nodes. We explicitly represent the nodes of the graph with a vector 
because rules may add and delete nodes, and thus we mark the existing nodes with a 1 
in the corresponding position of the vector. Although nodes and edges can be assigned 
a type (as in 1 14 1), here we omit it for simplicity. 

A production, or rule, p : L ^ Ris a partial injective function of simple digraphs. 
Using a static formulation, a rule is represented by two simple digraphs that encode the 
left and right hand sides. 

Definition 1 (Static Formulation of Production). A production p : L ^ R is stati- 
cally represented as p = (L = (L^ , L^); R = {R^ , R^)), where E stands for edges 
and V for vertices. 

A production adds and deletes nodes and edges; therefore, using a dynamic for- 
mulation, we can encode the rule's pre-condition (its LHS) together with matrices and 
vectors to represent the addition and deletion of edges and nodes. 

Definition 2 (Dynamic Formulation of Production). A production p : L ^ R is 
dynamically represented as p = {L = (L^ , L^); ,r^;e^ , r^), where and 
are the deletion Boolean matrix and vector, and are the addition Boolean matrix 
and vector (with a 1 in the position where the element is deleted or added respectively). 

The output of rule p is calculated by the Boolean formula R = p(L) = r v eL, 
which applies to nodes and edges (the a (and) symbol is usually omitted in formulae). 
Example. Fig. [T] shows an example rule and its associated matrix representation, in its 
static (right upper part) and dynamic (right lower part) formulations. ■ 

In MGGs, we may have to operate graphs of different sizes (i.e. matrices of different 
dimensions). An operation called completion 1 12 1 rearranges rows and columns (so that 
the elements that we want to identify match) and inserts zero rows and columns as 
needed. For example, if we need to operate with graphs Li and Ri in Fig.[T] completion 
adds a third row and column to R^ (filled with zeros) as well as a third element (a zero) 
to vector R^ . 

Compatibility. A graph (A/, V) is compatible if M and V define a simple digraph, 
i.e. if there are no dangling edges (edges incident to nodes that are not present in the 
graph). A rule is said to be compatible if its application to a simple digraph yields a sim- 
ple digraph (see ifTSl for the conditions). A sequence of productions s„ = p„; . . . ;pi 
(where the rule application order is from right to left) is compatible if the image of 
Sm = Pm', . . . ;pi is compatible, Vm ^ n. 
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Fig. 1. Simple Production Example (left). Matrix Representation, Static and Dynamic 
(right). 



Nihilation Matrix. In order to consider the elements in the host graph that disable a 
rule application, rules are extended with a new graph i^rQlts associated matrix specifies 
the two kinds of forbidden edges: those incident to nodes deleted by the rule and any 
edge added by the rule (which cannot be added twice, since we are dealing with simple 
digraphs) 

According to the theory developed in ifTSl . no extra effort is needed from the gram- 
mar designer to derive the nihilation matrix, ?& K = p (U) with Z) = ® , where 
(X) is the tensor product, which sums up the covariant and contravariant parts and mul- 
tiplies every element of the first vector by the whole second vector U4,l . Transposition 
will be represented by *. Please note that given an arbitrary LHS L, a valid nihilation 
matrix K should satisfy L^K = 0, that is, the LHS and the nihilation matrix should 
not have common edges. 

Example. The left of Fig. |2] shows, in the form of a graph, the nihilation matrix of the 
rule depicted in Fig. [T] It includes all edges incident to node 3 that were not exphcitly 
deleted and all edges added by pi. To its right we show the full formulation of pi which 
includes the nihilation matrix. ■ 
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Fig. 2. Nihilation Graph (left). Full Formulation of Production (center). Evolution of K 
(right). 



' In 1151 , K is written A''^ and Q is written Nr. We shall use subindices when dealing with 
sequences in Sec. [7] hence the change of notation. In the definition of production, L stands for 
left and R for right. The letters that preceed them in the alphabet (K and Q) have been chosen. 

^ Nodes are not considered because their addition does not generate conflicts of any kind. 



As proved in IfTSl (Prop. 7.4.5), the evolution of the nihilation matrix is fixed by the 
production. If R = p{L) = r veL then 

Q =p-\K) = evrK, (1) 

being Q the nihilation matrix of the right hand side of the production p. Hence, we have 
that (i?, Q) = {p{L),p^^{K)). Notice that Q ^ D in general though it is true that 

Example. The right of Fig. |2]shows the change in the nihilation matrix of pi when the 
rule is applied. As node 3 is deleted, no edge is allowed to stem from it. Self-loops from 
nodes 1 and 2 are deleted by p so they cannot appear in the resulting graph. ■ 

In lfT2l we introduced a functional notation for rules inspired by the Dirac or bra-kejl 
notation |[T]. Thus, we can depict a rule p : L ^ Ras R = p{L) = (L.p^, splitting the 
static part (initial state, L) from the dynamics (element addition and deletion, p). Using 
such formulation, the ket operators (i.e. those to the right side of the bra-ket) can be 
moved to the bra (left side) by using their adjoints. In this work we recast this notation 
more properly through MCL. 

Direct Derivation. A direct derivation consists on applying a rule p: L — > i? to a 
graph G, through a match m: L ^ G yielding a graph H. In MGGs we use injective 
matchings, so given p : L ^ R and a simple digraph G any m : L ^ G total injective 
morphism is a match for p in G. The match is one of the ways of completing L in G. 
In MGGs we do not only consider the elements that should be present in the host graph 
G (those in L) but also those that should not be (those in the nihilation matrix, K). 
Hence two morphisms are sought: . L ^ G and niK : K ^ G, where G is the 
complement of G, which in the simplest case is just its negation (see 1141151 ). 

Definition 3 (Direct Derivation). Given rule p : L ^ R and graph G = {G^ ,G^) 
as in Fig. \3][aj, d = (jp, m) — with m = {rriL, mx) - is called a direct derivation with 
result H = p* (G) if the following conditions are fulfilled: 

1. There exist rriL : L ^ G and ttik '■ K —>■ G^ total injective morphisms. 

2. mL{n) = rriKin), "in 6 . 

3. The match ttil induces a completion of L in G. Matrices e and r are then com- 
pleted in the same way to yield e* and r*. The output graph is calculated as 
H =p*{G) = r* V e*G. 

Remarks. The square in Fig. [3] (a) is a pushout. Item 2 is needed to ensure that L and 
K are matched to the same nodes in G. 

Example The right of Fig.[3]depicts a direct derivation example using rule pi shown in 
Fig. [T] which is applied to a graph G yielding graph H. A morphism from the nihila- 
tion matrix to the complement of G, rriK : K ^ G, must also exist for the rule to be 
applied. ■ 

Analysis Techniques. In II12I13I14I15I we developed some analysis techniques for 
MGGs. One of our goals was to analyze rule sequences independently of a host graph. 

We have followed the mathematical style instead of the one commonly used in physics, which 
should have been (p \ V). 
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Fig. 3. Direct Derivation (left). Example (right). 



For its analysis, we complete the sequence by identifying the nodes across rules which 
are assummed to be mapped to the same node in the host graph (and thus rearrange the 
matrices of the rules in the sequences accordingly). Once the sequence is completed, 
our notion of sequence coherence |12| allows to know if, for the given identification, 
the sequence is potentially applicable (i.e. if no rule disturbs the application of those 
following it). 

Given a completed sequence, the minimal initial digraph (MID) is the smallest graph 
that allows applying it. Conversely, the negative initial digraph (NID) contains all el- 
ements that should not be present in the host graph for the sequence to be applicable. 
Therefore, the NID is a graph that should be found in G for the sequence to be appli- 
cable (i.e. none of its edges can be found in G). If the sequence is not completed (i.e. 
no overlapping of rules is decided), we can give the set of all graphs able to fire such 
sequence or spoil its application. 

Other concepts aim at checking sequential independence (i.e. same result) between 
a sequence of rules and a permutation of it. G-congruence detects if two sequences 
(one permutation of the other) have the same MID and NID. It returns two matrices and 
two vectors, representing two graphs, which are the differences between the MIDs and 
NIDs of each sequence respectively. Thus if zero, the sequences have the same MID 
and NID. Two coherent and compatible completed sequences that are G-congruent are 
sequential independent. 

All these concepts have been characterized using operators A and y. They extend 
the structure of sequence, as explained in fl5l. Their definition is included here for 
future reference: 



tl / tl 



A*; {F{x,y))= V [/\{F{x,y))\ 

tl / y ^ 
^\l{G{x,y))= V /\{G{x,y)) 



(2) 



(3) 



y=ta \x=to 



Productions are the building blocks of sequences and sequences are the basic con- 
struction to study graph dynamics. All these concepts are further studied and general- 
ized in the present contribution. 

Some other important notions such as application conditions, graph constraints or 
reachability are just sketched or not even mentioned and left for further research. 



3 MCL: Monotone Complex Logic 



In this section we introduce Monotone Complex Logic (MCL), Preliminary Monotone 
Complex Algebra (PMC A) and Preliminary Monotone Matrix Algebra (PMMA). The 
term "logic" in the title should be understood as in fist-order logic or propositional 
logic (arguably, "calculus" might be more appropriate). It has been called complex to 
resemble the similarities with complex numbers and how they are defined out of the 
real numbers. Monotone because we are not defining the negation of complex terms 
(see below). 

Monotone complex logic is in our opinion of interest by itself, but it is introduced 
here due to its usefulness for Matrix Graph Grammars (MGGs). First, it permits a com- 
pact reformulation of grammar rules. Second, the numerical representation that will be 
introduced in Def.[8]establishes a link between graphs in MGGs and Q[«], although the 
operations we are interested in are not addition and multiplication. Also, any production 
p induces the evolution of a pair of graphs {L.K) {R,Q) = {p{L),p-\K)). Pro- 
ductions will be reinterpreted by encoding them as complex formulas and representing 
their actions as a "Hermite product". MCL will allow us to measure the size of graphs 
via a natural norm. Finally, sequential notions of MGGs such as independence, initial 
digraphs, coherence, etcetera, will be thus recasted and extended. 

Definition 4 (Complex Formula). A complex formula z = {a,b) consists of a cer- 
tainty part 'a' plus a nihil part 'b', where a and b are propositional logic formulas us- 
ing adjacency matrices as propositional variables. Two complex formulas zi = (ai,6i) 
and Z2 = (02, 62) are equal, zi = Z2, if and only if ai = a2 and bi = 62. 

Monotone Complex Logic is the formal system whose propositional variables are 
complex formulas with logical connectives { v , a }. We will not go further because we 
are more interested in an algebraic development of the theory. 

Throughout the present contribution, complex formula, complex term and Boolean 
complex will be used as synonyms. Next, some basic operations on Boolean complexes 
are introduced. 

Definitions (Basic Complex Operations). Lef z = (a, 5), zi = (ai,6i) and Z2 = 
(02, ^2) be complex terms. The following operations are defined componentwise: 

- Addition: zi v 22 = (o-i v 02, 5i v 62)- 

- Multiplication: z\ /\ Z2 = z\ Z2 = (ai 02 v 61 62, Qi ^2 v 0261). 

- Conjugation: z* = {b,a). 

- Dot Product: (zi, 22) = zi z*. 

The notation •) is used for two reasons. First, we would like to highlight the 
similarities with scalar products. There is however no underlying linear space so this 
is just a convenient notation. Second, we will see that it coincides with the functional 
notation introduced in 1 12 15 1. 

The dot product of two Boolean complexes is zero (orthogonal) if and only if any 
element of the first complex term is included in both the certainty and nihil parts of 
the second complex term. Otherwise stated, if zi = (ai,6i) and Z2 = (02,62), then 



(zi, z-i) = <f» a\a2 = 0162 = ^1^2 = ^1^2 = 0. Let's say that a < b if ab = a, i.e. 
whenever a has alb also has a 1 (graph a is contained in graph 6). Previous identities 
can be rephrased as ai < 02, ai < 62, 61 < 02 and 5i < 62. This is equivalent to 
(oi V 61) < (0262). Orthogonality is directly related to the common elements of the 
certainty and nihil parts. 

A particular relevant case is when we consider the dot product of one element z = 
(a, b) with itself. In this case we get (a v b) < (ab), which is possible if and only if 
a = b. We shall come back to this issue. 

Definition 6 (Preliminary Monotone Complex Algebra, PMCA).r/ie5efG5' = {z\z 

is a complex formula] together with the basic operations introduced in Def. \5\will be 
known as preliminary monotone complex algebra. 

We will get rid of the term "preliminary" in Def. [TTl when not only the adjacency 
matrix is considered but also the vector of nodes that make up a simple digraph. 

We introduce a subalgebra of the preliminary monotone complex algebra to be 
known as preliminary monotone matrix algebra (PMMA). It is useful due to its re- 
lationship with MGGs. 

Definition? (Preliminary Monotone Matrix Algebra, PMMA). Lef = (ai,bi) e 
&. Define the equivalence relation zi ~ Z2 3c = (c, c) | a2 = aiw candb2 = 61 vc. 
Then, 

io' = ©7 ^ (4) 
is the preliminary monotone matrix algebra. 

A graph and any of its possible nihilation matrices in MGG do not share any edge 
(L^K = 0, refer to [15]). So when representing the left hand side of a production 
in MCL, its complex term is made of a digraph in the certainty part and some valid 
nihilation matrix in the nihil part. Intuitively, PMMA is made of the valid complex 
terms in this sense (i.e. those that do not share any edge). 

It is not difficult to check reflexivity, symmetry and transitivity for ~. The equiva- 
lence relation permits the simplification of those elements that appear in both the cer- 
tainty and nihil parts (eliminating non- valid complex terms). 

A more complex-analytical representation can be handy in some situations and in 
fact wiU be preferred for the rest of the present contribution: z = (a, b) 1 — > z = avib. 
Define one element i, that we will name nil term or nihil term, with the property i Ai = 
1, being i itself not equal to 1. Then, the basic operations of Def.|5l following the same 
notation, can be rewritten: zi v Z2 = (ai v a2) v i{bi v 62), Z1Z2 = (ai v ibi) a 
(02 V ib2), z* =b V ia and the same for the dot product. 

Notice that the conjugate of a complex term z e <&' that consists of certainty part 
only is z* = (a v iO)* = 1 via. Similarly for one that consists of nihil part alone: 
z* = {0 V ib)* = b V i. If z e i)' then they further reduce to a v iO and v ib, 
respectively, i.e. they are invariant^ Also, the multiplication reduces to the standard 
and operation if there are no nihil parts: (ai v iO){a2 v iO) = aia2- 

* Notice that 1 v ia = (a v a) v ia = a v iO and 5 v il = fo v i (6 v 6) = v i6. 



Proposition 1. Let x,y,z e ©' and z\,Z2 6 fi' . Then, (x v y,zy = (x, z) v (y, z), 

<zi, Z2> = <Z2, zi>* anc/ (2:12:2)* = -2* 2|. 

f roof 

□The first identity is fulfilled by any complex term and follows directly from the defi- 
nition. The other two need the equivalence relation (simplification), i.e. they hold in fy' 
but not necessarily in For the second equation just write down the definition of each 
side of the identity: 

<zi, 22> = {aj)2 V 02^1) V i {aia2 v 6162) 

<Z2,2i>* = [aife va26i v (ai6i va2 62)] v i [aia2 v 6162 v (ai6i va2 52)]. 

Terms ai6i v 0,2^2 vanish as they appear in both the certainty and nihil parts. The 
third identity is proved similarly. ■ 

Notice however that (21 v 22)* z^wz*. It can be checked easily as (21 v 22)* = 
[(fli V 02) V i (61 V 62)]* = 6162 vi aia2 but zf v z^ = (61 v 62) vi (oi v 02). This 
impUes that, although (21 v 22, 2) = (21, 2) v {22, 2), we no longer have sesquilineal- 
ity, i.e. it is not linear in its second component taking into account conjugacy: 

2 [(61 V 62) V i (ai V 02)] = <2, 2l V 22> # <2, 2l> V {z, 22> = 2 [6162 V 10.10,2^ ■ 

In fact the equality <2, 21 v 22) = {z, 21) v <2, 22) takes place if and only if 21 = 22. 

4 Numerical Representation, Norm and Distance 

This section introduces an application £ that assigns a complex number in the unit inter- 
val C ([0, 1]) to any Boolean complex. It is not a homomorphism: neither £{x v y) = 
£{x) + £{y) nor £{xy) = £{x)t{y) hold. Application £ provides some geometric in- 
tuition. A norm and a conditional norm are defined out of the dot product of Def. |5j 
Finally we will define the distance between two complex terms. 

Definition 8 (Rational Enconding). Let g = (gj) . ^^^^ be a simple digraph. Its 
rational encoding is given by 

£{g) = 2 (2-'^g^ , (5) 

k = l 

where ik = \^] and jk = k - n[^\. 

It has become customary to represent the lowest integer above m as [m] and the 
biggest integer below m as [mj. These functions are known as ceiling and floor, respec- 
tively. The indices ik and jk are just the integer quotient and the remainder. They are a 
convenient way to visit all the elements of the adjacency matrix ordered by columns. 

As the elements of &' are adjacency matrices we can define G' = £ {&)■ If 2 = 
(a, b) 6 25', then £{z) = £ia) + i£{b). Analogously, W = £ {S)'). 



Example. £(Li) = 0'10101l2 = 0'671875io, where Li is the LHS of pi in Fig.[I] If 
then £{£) = O'Ollz + iO'h = 0'375io + iO'Sio- Subindices 

indicate the numbering system: base 2 or base 10. ■ 

The rational encoding is similar to the standard one given in the literature (which we 
will call natural encoding). The main difference is that we use negative powers. This, 
in the limit, would make us consider the interval [0, 1] as the underlying space instead 
of the natural numbers N. 

In the present contribution we will deal only with finite graphs. Hence, the codomain 
of £ is [0, l[nQ, in fact only terminal rationalslf] Nonetheless, a small digression seems 
appropriate as at some point in the future we will be interested in the asymptotic be- 
haviour of algorithms and, hence, it will be more convenient to consider the interval 
[0, 1]. Note that is the symmetric element with respect to the v operation. This is the 
graph with no edges. It should be also natural to include 1 because it should be the neu- 
tral or identity element with respect to the a operation (the graph that has any possible 
edge) if the assumption 1 = 0'T2 was made. However, this assumption is not adequate 
for MGGs. For example, if we want to consider the graph g that only has the self-edge 
(1,1) and no other one, and the number of nodes of the graph was countable, we would 
be tempted to write g = 0'l2 v 0'0T2. But 0'0T2 = 0'l2 with the standard operations in 
C. We would thus be asking for the self-edge (1, 1) to be present and not to be present, 
according to the interpretation of the certainty and nihil parts in MGG. Seemingly, the 
Archimedean property fails as for graphs there are nontrivial infinitesimals. This im- 
pUes for example that any diagonal in Fig.|4]does not belong to the MGG characteristic 
set. 

Throughout the present contribution we may make the following abuse of notation. 
Let A be the adjacency matrix of some simple digraph. When we make the operation 
A V A we do not obtain 1 but 1a, which is the matrix with all its elements 1 according 
to vector , i.e. Aw A = AX (x)A^. This has a clear relationship with the definition of 
the nihilation matrix in ifTSl . See the definition of conditional norm below. Also, we can 
interpret \a as the characteristic or indicator function of all edges potentially incident 
to nodes of digraph A (the smallest complete digraph that contains A). The notation 
Xa{x) is also standard. It is defined xa{x) = 1 if x e A^ ® A^ and zero otherwise. 

The rational encoding in Def.[8]establishes an injection from the preliminary mono- 
tone complex algebra & into the complex numbersjj where the certainty part becomes 
the real part, the nihil part becomes the imaginary part and the nihil term becomes 
i = V— T. Figure |4]represents the preliminary monotone matrix algebra Sj' as a proper 
subset of C which is a well known fractal. 



^ By hypothesis, binary representation of rationals will always be terminating because any 
dyadic rational number 1/2" has a terminating binary numeral (although other rational num- 
bers recur). To fix the notation: a terminal number is O'OllOl while a recurring number is 
o'0n = O'OllOllOll . . . 

^ Though, as commented in this section, just as sets. The morphism is an injection if the number 
of nodes is fixed. For example, 0'10102 may be the graph with two nodes and edges (1, 1) and 
(1, 2) or 0'1010000002, the graph with three nodes and edges (1, 1) and (3, 1). 
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Fig. 4. Rational Encoding of the Preliminary Matrix Algebra as a Subset of C 



Proposition 2. The characteristic function of the rational encoding of the preliminary 
matrix algebra approaches the Sierpinski gasket as the number of nodes increases. 

Proof 

□The set corresponds to the zeros of the and function (coloring zeros in black). It is the 
Sierpinski gasket due to the Lucas correspondence theorem |5| (see also L18|) which 
can be used to compute the binomial coefficient (^) mod 2 with bitwise operations: 
LaK. This tells us that the parity of the function (^) (this is what the function mod 2 
does) is the same as that of L aK. In our case L is the abscissa. Its negation just reverts 
the order (it is a symmetry) and does not change the shape of the figure. As commented 
above in this section, we do not want the diagonals to belong to the set in the limit. ■ 

The dot product introduced in Sec.[3]induces a "norm" || • || in ©'. It is not a norm in 
the sense of linear algebra (no underlying vector space) but shares some of its properties 
and will play a similar role: it will be used as a means to measure the size of a graph. 



Definition 9 (Norm - Conditional Norm). Let y,z e &. Its norm is the application 
IHI : & ->■ [0, 1] definedby 

\\z\\=e{(z,zy). (6) 

The conditional norm of z with respect to y is given by 

ll^ll - My - (7) 

" \\y\\ 

The following identities show that the norm (before applying £) returns what one 
would expect. Equation ( fTOl i is particularly relevant as it states that in MCL the cer- 
tainty and nihil parts are in some sense mutually exclusive, which together with eq. 
(fTTI) suggest the definition of Sj' as introduced in Sec. [3] Notice that this fits perfectly 



well with the interpretation of L and K given in [ 15 1. 

<(a,0),(a,0)>= (a,0)(l,a) = (a 1 v a, a a v 1) = a (8) 
<(0, 6), (0, 5)> = (0, 6) (5, 1) = (5 1 V 5, 56 V 1) = 5 (9) 
{(c, c), (c, c)) = (c, c) (c, c) = (cc V CC, CC V cc) = 0. (10) 

The dot product of one element with itself gives rise to the following useful identity: 

(z, z) = z z* = {ah V a5) v i (bb v aaj = a@b. (11) 

Equation (fTTT l admits two readings. In first place, it tells how to factorize one of the 
basic Boolean operations, xor. Secondly, and more relevant to us, it justifies the use of 
xor as a norm for complex termsQ Besides, it follows directly from the definition that 
||z||=land||z*|| = ||z||. 

The conditional norm || • | y 1| reduces the total space to graph y.lf z = ai vibi e & 
and y = a2 V ib2 e after some simple algebraic operations 

t (02 ©52) 

is obtained. If we further have that y, z e Sj', then 

. ^((^i v ^i)(«2 v52)) ^ ^^3^ 

I (02 V 62) 

In particular this implies that ||2:|?;|| = 1 <f» y is a subgraph of z. The following propo- 
sition highlights the similarities between norms in linear spaces and the one introduced 
inDef.|9] 

Proposition 3. Let y, z e 9j' and zi, Z2 e & . Then, for norms and conditional norms 
we have that 

||z||=0^z = 0, (14) 

lly^ll = lly|| A ||z||, (15) 

Ikl VZ2IK pill V ||Z21|. (16) 

Proof 

□Identity ( fT4b is derived from ( ITTT i. Some simple manipulations prove (fTsT i (where and 
is performed bitwise). Inequality (fTST i is not difficult either where, again, or on the right 
hand side is applied bitwise: 

V IIZ2II = i{aj)i) V £{a2b2) v -^(0151) v £(0262) 
\\zi V 2211 = ^(ai5i52) V i{a2bib2) v £{biaia2) v ^(620102). 

Comparing term by term, the inequality follows. For example, there must be at least the 
same numbers of I's in ai5i than in ai5i62, so £{aibi) ^ ^(ai5i52).B 



' Recall that in complex analysis we have z z* 



The xor operation z = zi @ Z2 is the number of distinct elements in zi and Z2. 
The number of ones that appear in z tells the number of atomic operations that must be 
performed in order to transform zi in Z2 (or viceversa) in the sense of MGG produc- 
tions. Therefore, it seems natural to define the distance between two complex terms as 
follows: 

Definition 10. Let Zj = aj v ibj with Zj e &, j = 1,2, and define wi = ai v ia2 and 
W2 = bi V ib2- Then, the distance between zi and Z2 is given by 

d{zi,Z2) = \\WI\\@\\W2\\. (17) 

It is an easy exercise to check that d fulfills the axioms of a metric: d{x, y) ^ 0, 
d{x, y) = X = y, d{x, y) = d{y, x) and the triangle inequality d{x, z) ^ d(x^ y) + 
d{y, z). The triangle inequality follows from d{x, z) = d{x, y) @ d{y, z) and the fact 
that Va > 0, V6 > we have that a @ 6 ^ a + 6. See for an application of the xor 
metric. 

Notice that in 9)' only one of the terms contributes to the norm if Z2 = pizi) for 
some production because the actions on the certainty part completely determine the 
nihil part, as proved in [15J (Prop. 7.4.5). 

5 Production Encoding 

In this section we introduce the Monotone Complex Algebra and the Monotone Matrix 
Algebra, that not only consider edges but also nodes. Compatibility issues may appear 
so we study compatibility for a simple digraph and also for a single production (com- 
patibility for sequences will be addressed in Sec. |7]i. Next we turn to one of the main 
topics in this paper: how to characterize MGG productions using MCL and the dot 
product of Def. |5] The section ends introducing swaps and providing some geometric 
interpretations. 

To get rid of the "preliminary" term in the definitions of ©' and S)' (Defs.|6]and|7l 
resp.) we shall consider an element as being composed of a matrix term and a vector 
of nodes. Hence, we have that £ = (L^ v iK^, v iK^^ where E stands for edge 
and V for verfex0 Notice that v iK^ are matrices and v iK^ are vectors. 

Definition 11 (Monotone Complex and Matrix Algebras). © and (the Monotone 
Complex and Matrix Algebras, resp.) are defined as their preliminary counterparts - 
see Defs.^and^- but considering elements of the form: 

C= {L^ w iK'^, V iK^) . (18) 

We also introduce G = £ (6) and M = £ (S)). 

Concerning ©, a production p : (3 ^ (& consists of two independent productions 
P = (Pi)P2) - being pi MGG productions as those introduced in ifTSll - one acting on 
the certainty part and the other on the nihil part: 

n = p{C)=pc{L)wipN{K)=RwiQ. (19) 

* If an equation is applied to both edges and nodes then the superindices will be omitted. They 
will also be omitted if it is clear from context which one we refer to. 



As there are no restrictions on pc and pat if we stick to 25, it is true that V(?i , 32 6 ©, 3p 
such that p{gi) = 32- Recall from Sec. |2]that productions in MGG have as domain 
and codomain. Moreover, they must fulfill pjy = p^ . Unless otherwise stated, we will 
concentrate on MGG productions for the rest of the paper. 

We want PAT to be a production so we must split it into two parts: the one that acts 
on edges and the one that acts on vertices. Otherwise there would probably be dangling 
edges in the nihil part as soon as the production acts on nodes. The point is that the 
image of the nihil part with the operations specified by productions are not graphs in 
general, unless we restrict to edges and keep nodes apart. This behaviour is unimportant 
and should not be misleading. 
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Fig. 5. Potential Dangling Edges in the Nihilation Part 



Example. nTo the left of Fig. |5] we have drawn the certainty part of a production p 
that deletes node 1 (along with two incident edges) and adds node 3 (and two incident 
edges). Its nihil counterpart for edges is depicted to the right of the same figure. Notice 
that node 1 should not be included in K because it appears in L and we would be 
simultaneously demanding its presence and its absence. Therefore, edges (1, 3), (2, 1) 
and (3, 1) - those with a red dotted line - would be dangling in K (red dotted edges do 
belong to the graphs they appear on). The same reasoning shows that something similar 
happens in Q but this time with edges (1, 3), (3, 1), (3, 2) and (3, 3) and node 3. 

This is the reason to consider nodes and edges independently in the nihil parts of 
graphs and productions. In K, as nodes 1 and 3 belong to L, it should not make much 
sense to include them in K too, for if K dealt with nodes we would be demanding their 
presence and their abscense. In Q the production adds node 3 and something similar 
happens. ■ 

Now that nodes are considered compatibility issues in the certainty part may show 
up. The determination of compatibility for a simple digraph under MCL is almost 
straightforward. Let g = (g^ v ig^,gc v ig^) G -Q. Potential dangling edges are 

given by Dg = g^ ® g^, so the graph g will be compatible if g^ Dg = 0. If g e 53 
there are no common elements between the certainty and nihil parts and Dg < gfj. 

A production p{C) = p{L v iK) = R v iQ = TZ is compatible if it preserves 
compatibility, i.e. if it transforms a compatible digraph into a compatible digraph. This 
amounts to saying that RQ = 0. 

Recall from Sec. |2] that grammar rules actions are specified through erasing and 
addition matrices, e and r respectively. Because e acts on elements that must be present 



and r on those that should not exist, it seems natural to encode a production as 

p = e V ir. (20) 

Our next objective is to use the dot product - see Def.|5]- to represent the application 
of a production. This way, a unified approach would be obtained. To this end define the 
operator P : ^ © by 

p = e w ir I — > P{p) = er w i {e w r) . (21) 

Proposition 4 (Production). Let C and TZ be the left and right hand sides, resp., as in 
Def. [77] eq. Iil8\l . and P as defined in eq. (121b . Then, 

n = {C,P{p)). (22) 

Proof 

□The proof is a short exercise that makes use of some identities which are detailed right 
afterwards: 

<AP(p)> = <(i,i^),(er,e vr)> = 

= (erL V (e V r)K,erK v (e v r)L) = 

= {r weL.ewrK) = {p{L),p-\K)) = TZ. (23) 

Apart from equation (4.13) of Prop. 4.1.4 in ifTsl which states that rL = L, we have 
used the following identities: 

(e V r)K = eK v rK = rK = r(r v eD) = r. 

erK = r(er v eeZ3) =rK. 
(e V r)L = eL w rL = eL = e. 

We have also used that re = r (again Prop. 4.1.4 in jTSll ). rD = r due to compatibility 
and rL = almost by definition. Besides, Prop. 7.4.5 in ifTSI has also been used, which 
proves that Q = p^^ (^)-' 

The production is defined through operator P instead of directly asp = ervi{evr) 
for several reasons. First, eq. ( |20| ) and its interpretation seem more natural. Second, 
P{p) is self-adjoint, i.e. P(p)* = P(p), which in particular implies that ||P(p)|z|| = 1, 
Vp, being z the "total" graph, i.e. the graph with respect to which completions are per- 
formed (for completion refer to lITSI . Sec. 4.2). Therefore, the norm would not measure 
the size of productions (interpreted as graphs according to eq. ( |20] |) and we would be 
forced to introduce a new norm. This is because 

(P(p), P(p)) = (er V i{e v r)) (er v i{e v r))* = er v e v r = 1^. 

By way of contrast, |jp|j = e@r = evr. With operator P the size of a production is the 
number of changes it specifies, which is appropriate for MGGs0 Moreover, due to eq. 



' One of our objectives is to look for an appropriate measure of the number of actions that would 
eventually transform one graph into another. 



( |20l ). geometrically, grammar rules are also the Sierpinski gasket. In fact, the codomain 
of operator P is the diagonal e + r = 1 (refer to Figs.|4]and|7]l. 

Complex logic encoding puts into a single expression the application of a grammar 
rule, both L and K. Also, it links the functional notation introduced in 1.15,1 and the dot 
product of Sec. [3] 

Theorem 1 (Surjective Morphism). There exists a surjective morphism from the set 
ofMGG productions on to the set of self- adjoint graphs in S^. 

Proof 

□It is not difficult to check that z is self-adjoint if and only if 1 1 z 1 1 = 1 z : on the one hand, 
if z = a V ia then (z, z) = zz* = (a v ia){a v ia) = ava = Iz- On the other hand, if 
we have z = a v ib and ||z|| = a®b = Iz then a = b. Note that ||_zl| = 1^ is equivalent 
to asking for the conditional norm to be equal to 1 with respect t(0 Z = z®z^: 
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The surjective morphism is given by operator P. Clearly, P is well-defined for any 
production. To see that it is surjective, fix some graph g = gi v ig2 such that \\g\\ = Ig. 
Then, g = gi v ig^. Any partition ofgi as or of two disjoint digraphs would do. Recall 
that productions (as graphs) have the property that their certainty and nihil parts must 
be disjoint. ■ 

The operator P is surjective but not necessarily injective. It defines an equivalence 
relation and the corresponding quotient space. We will be led to a reinterpretation of 
the notion of production in Matrix Graph Grammars. 

Definition 12 (Swap). The swap space is defined as W = ?j/P(Sj). An equivalence 
class in the swap space will be called a swap. The swap w associated to production 



p : Sj ^ is w 



P{p), i.e. p € Sj 
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Fig. 6. Example of Productions 



Example. nLet p2 and be two productions as those depicted in Fig.|6] Their images 
in H are: 



P(P2) = P{p3) 



"0 0" 




"1 1" 
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V i 
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= w 



(24) 



The tensor (Kronecker) product in this contribution will always be used on nodes. The V 
superindex will be omitted in this case: z^z^ = 2^ (g) (z*)^, where * stands for transposition. 
" Acording to eq. J20t . any element in S) can be interpreted as a production and viceversa. 



They appear to be very different if we look at their defining matrices L2 , and 
i?2 , i?3 or at their graph representation. Also, they seem to differ if we look at their 
erasing and addition matrices: 





"10" 




"11" 




"0 1" 




"00" 


62 = 


1 


63 = 


00 


f 2 = 


00 


ri, = 


1 



However, they are the same swap as eq. (l24l i shows, i.e. they belong to the same 
equivalence class. Notice that both productions act on edges (1,1), (2,2) and (1,2) 
and none of them touches edge (2, 1). This is precisely what eq. (l24l l says as we will 
promptly see. 

Swaps will be of help in studying and classifying the productions of a grammar. 
For example, there are 16 different simple digraphs with 2 nodes. Hence, there are 256 
different productions that can be defined. However, there are only 16 different swaps. 
From the point of view of the number of edges that can be modified, there is 1 swap 
that does not act on any element (which includes 16 productions), 4 swaps that act on 1 
element, 6 swaps that act on 2 elements, 4 swaps that act on 3 elements and 1 swap that 
acts on all elements. ■ 

There is a simple geometrical interpretation of P in terms of the rational encoding 
of the productions p and the way they are transformed. Let the principal diagonal be the 
closest lin^to x + y = 1. Operator P assigns the same element P{p) in the principal 
diagonal to any element p of a parallel line (to the left of Fig. [7] three parallel diagonals 
are represented ). P{p) has the same certainty part as the biggest certainty part of any 
p in the parallel line. See Fig. |7] for the transformation of three sets of productions 
(circles) into their associated swaps (squares). Geometrically, P can be thought of as 
the composition of two projections: one along the corresponding diagonal and another 
parallel to the abscissa axis. 




Fig. 7. Rational Encoding of the Transformation of Productions via Operator P 



"Closest" because the line x + y = 1 is approached as the number of nodes in the graphs tend 
to 00. 



The name swap has been chosen because of the way they act on elements in f). 
Following previous example, let's consider the swap w = O'OlOO v iO'lOll. Let's 
also consider a generic element C = Q'L\L\L\L^ v iO'KlKfK2K2- The image 
w(£) = Q' KIL\K2K2 V iO'L\KfLlLl swaps the elements that appeal" in the nihil 
part of w and keeps unaltered those that appear in its certainty part. Swaps summarize 
the dynamics of a production, independently of its left hand side. Notice that, because 
swaps are self-adjoint, it is enough to keep track of the certainty or nihil parts. So one 
production is fully specified by, for example, its left hand side and the nihil part of its 
associated swapPi 

We can reinterpret actions specified by productions in Matrix Graph Grammars un- 
der MCL: instead of adding and deleting elements, they interchange elements between 
the certainty and nihil parts. 

The geometrical intepretation of the actions of swap tui = P{p) = 0'0T2 v ?0'l2 
(first element is swapped and the rest remain unaltered^ is a reflection with respect to 
the diagonal y = x in the Sierpinski gasket S. See the right of Fig. [T] We are consid- 
ering those Boolean complexes that are complementary in the places where the swap 
interchanges elements^ Swap W2 = 0'10T2 v j0'0l2 divides S into two regions, 5*1 
and S2, along the line y = x. It acts again as a reflection, but independently in Si and 
5*2. The reflection is defined with respect to the line y = x — ^ in Si and with respect 
to y = X + i in 5*2. The swap W12 = 0'00T2 v iO'Ol is the composition of wi and W2. 
The order does not matter. Geometrically, corresponds to a reflection with respect 
to y = X and another reflection with respect to y = x + ^ (only one of them applies). 




Fig. 8. Image of Simple Digraphs via Swaps 



Given a swap and a complex term jC, it is straightforward to calculate the production having L 
as left hand side and whose actions agree with those of the swap. 

The notation T stands for "as many ones as nodes". In this contribution we deal with finite 
graphs. 

For example, if wi acts on zi = O'OllO v iO'OOOl then toi(zi) = zi.lf the first element in 
zi is the same in both graphs (in fact zero) then zi is a fixed element of wi. 



Let's consider those digraphs L whose certainty and nihil parts are complementary. 
Figure |8] has L on the x-axis and the certainty part of the swaps w in the y-axis. The 
z-axis represents its image w{L) = w @ L. Both figures are the same surface, but the 
second has been rotated 7r/4 radians clockwise around the z-axis. Notice that the more 
elements of L that we fix as zero (instead of as the complement of its certainty part) the 
closer w @L gets iow /\ L. 

6 Coherence and Initial Digraph 

So far we have extended MGGs by defining the transformations in (5 and The the- 
ory will be more interesting if we are able to develop the necessary concepts to deal 
with sequences of applications rather than productions alone. Among the two most ba- 
sic notions are coherence and the initial digraph. We shall reformulate and extend the 
concepts introduced in \ \5\. 

Recall that coherence of the sequence s = p„; . . . ; pi guarantees that the actions of 
one production Pi do not prevent the actions of those sequentially behind it: pi+i , ■ ■ ■ ,Pn 
(the first production to be applied in s is pi and the last one is p„; the order is as in com- 
position, from right to left). 

Theorem 2 (Coherence). The sequence of productions s = pn', ■ ■ ■ pi is coherent if 
the Boolean complex C = C"*" v iC^ = 0, where 

n 

C+ = Y (i?, v"+i (e.ry) V Aj-^ (e, r,)) (25) 

and 

n 

C-=\/ (Qj V"+i (e, r,) V K, A{-' (vy e,)) . (26) 

Proof 

□C+ V iC^ = C+ = = 0. The certainty part C+ is addressed in fB] 
and C~ = can be proved similarly. The reader is invited to consult the proof of Th. 
4.3.5 in ifTsll plus Lemma 4.3.3 and the explanations that follow Def. 4.3.2 in the same 
reference. Next, a sequence of two productions s = P2;Pi is considered to show the 
way to proceed for . 

In order to decide whether the application of pi does not exclude p2 (regarding 
elements that appear in the nihil parts) the following conditions must be demanded: 

1 . No common element is deleted by both productions: 

6162 = 0. (27) 

2. Production p2 does not delete any element that the production pi demands not to 
be present and that besides is not added by pi: 



e2Kiri = 0. 



(28) 



3. The first production does not add any element that is demanded not to exist by the 
second production: 

riK2 = 0. (29) 



Altogether we can write 6162 v rie2Ki v riK2 = 62(61 v riKi) v riK2 = 
62Q1 V riK2 = 0, which is equivalent to 

627^2 Qi V eiriK2 = (30) 

due to basic properties of MGG productions (see Prop. 4.1.4 in ifTSl ). For a sequence 
that consists of three productions, s = P3;p2',Pi, the procedure is to apply the same 
reasoning to subsequences p2',Pi (restrictions on p2 actions due to pi) and P3;p2 (re- 
strictions on p3 actions due to pi) and or them. Finally, we have to deduce what has to 
be imposed on p^ actions due to pi, but this time taking into account that p2 is applied 
in between. Altogether: 

Qi (62 V ^263) V (5263 V K2ri V K3 (rie2 v 7-2) . (31) 

We may proceed similarly for four productions. Equation ( |26] ) can be deduced ap- 
plying induction on the number of productions. 

To see that eq. ( |30] | implies coherence we only need to enumerate all possible ac- 
tions on the nihil parts. It might be easier if we think in terms of the negation of a 
potential host graph to which both productions would be applied (G) and check that 
any problematic situation is ruled out. See table [T] where D is deletion of one element 
from G (i.e., the element is added to G), A is addition to G and P is preservation^ For 
example, action A2 ; A\ tells that in first place p\ adds one element e to G. To do so this 
element has to be in 61 (or be incident to a node that is going to be deleted). After that, 
P2 adds the same element, deriving a conflict between the rules. 
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Table 1. Possible Actions for Two Productions 



This proves G^ = for the case n = 2. When the sequence has three productions, 
s = P3',P2',Pi, there are 27 possible combinations of actions. However, some of them 
are considered in the subsequences P2', Pi and P3\P2- Table |2] summarizes them. 

There are four forbidden actionsO D3; Di, A3; Pi, P3; Di and A3] Ai. Let's con- 
sider the first one, which corresponds to rirs (the first production adds the element - it 
is erased from G - and the same forp3). In Table|2]we see that related conditions appear 

'* Preservation means that the element is demanded to be in G because it is demanded not to 
exist by the production (it appears in A'l) and it remains as non-existent after the application 
of the production (it appears also in Qi). 

" Those actions appearing in table [T] updated forps. 
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Table 2. Possible Actions for Three Productions 



in positions (1,1), (4, 1) and (7, 1). The first two are ruled out by conflicts detected in 
P2',Pi and P3;p2, respectively. We are left with the third case which is in fact allowed. 
The condition r^ri taking into account the presence of p2 in the middle in eq. (l3Tl i is 
contained in i^3rie2, which includes ^162^3. This must be zero, i.e. it is not possible 
for pi and p^ to remove from G one element if it is not added to G by p2- The other 
three forbidden actions can be checked similarly. 

The proof can be finished by induction on the number of productions. The induction 
hypothesis leaves again four cases: D„;Di, A„;Pi, Pn]Di and An;Ai. The corre- 
sponding table changes but it is not difficult to fill in the details. ■ 

There are some duplicated conditions, so it could be possible to "optimize" C. The 
form considered in Th. |2] is preferred because we may use A and y to synthesize 
the expressions. Notice that eq. (l27T i is already in C through eq. (IZST l. which demands 
eiL2 = (as 62 a L2 we have that eiL2 = ^ £162 = 0). Condition ( l28T l is 
e2Kifi = e2riri v e2rieiDi = 6261 Z?i, where we have used that Ki = p (Di)- 
Note that those eiDi ^ are the dangling edges not deleted by pi. Finally, equality 
(|29] | is riK2 = rip2 (D2) = fi {r2 v 6202) = rir2 v rie2l?2- The first term (rir2) 
is already included in C and the second term is again related to dangling edges. Poten- 
tial dangling edges appear in coherence and this may seem to indicate a possible link 
between coherence and compatibility0 
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Fig. 9. Example of Coherence 



Example. Let's consider the sequence s = ps; P4. Recall that the order of applica- 
tion is from right to left so p4 is applied first and p^ right afterwards. Let p^ and p^ be 



Compatibility for sequences is characterized in Sec.|7] Coherence takes into account dangling 
edges, but only those that appear in the "actions" of the productions (in matrices e and r). 



those productions depicted in Fig.|9] Once simplified, its coherence term is 



C(s) = C~^{s) V iC (s) = (-R4r5 v ^564) v i{Qie^ v K^i^r^) = 
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Coherence problems appear in this example for several reasons. Edge (2, 3) is added 
twice while self-loop (3, 3) is first deleted in p4 and then used in ^5. Edge (1, 3) be- 
comes dangling because production 735 deletes node 1. Edge (2,3) appears in C~{s) 
for the same reason that makes it appear in («)•■ 

The minimal initial digraph M(s) for a completed sequence s = p„; . . . ;pi was 
introduced in |15 | as a simple digraph that permits all operations of s and that does 
not contain a proper subgraph with the same property. The negative initial digraph has 
a similar definition but for the nihil part. Theorem |5] encodes as a complex term the 
minimal and negative initial digraphs, renaming it to initial digraph. 

Now we are interested in what elements will be forbid den an d which ones will be 
available once every production is applied^l Matrix D = e ® e* specifies what edges 
can not be present because at least one of their incident nodes have been deleted. Let's 
introduce the dual concept: 

T = a (e(x)e*) . (32) 

T are the newly available edges after the application of the production because of 
the addition of nodes0The first term, r ® r*, has a one in all edges incident to a vertex 
that is added by the production. We have to remove those edges that are incident to 
some node deleted by the production, which is what e®e} does. 



0011 0001 

0012 T=0012 
1113 0113 



Fig. 10. Available and Unavailable Edges After the Application of a Production 



Recall that whenever the tensor (Kronecker) product is used, we refer to the vector of nodes 
so the V superscript is omitted. For example R R' = R^ (g) (i?^)*. The * stands for 
transposition. 

This is why T does not appear in the calculation of the coherence of a sequence: coherence 
takes care of real actions (e, r) and not of potential elements that may or may not be available 
(D,T). 
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Example. oFigurelTOldepicts to the left a production q that deletes node 1 and adds node 
3. Its nihil term and its image are 
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To the right of Fig. [TO] matrix T is included. It specifies those elements that are not 
forbidden once production q has been applied. ■ 

Matrices D and T do not tell actions of the production to be performed in the com- 
plement of the host graph, G. Actions of productions are specified exclusively by ma- 
trices e and r. 



Theorem 3 (Initial Digraph). The initial digraph M{s) for the completed coherent 
sequence of productions s = pn; . . . ; pi is given by 

M (s) = Vi (rxLy V i ejT^Ky) . (33) 



Proof (sketch) 

□The proof proceeds along the lines of that for Th. 4.4.2 in 1 15 1, which in essence starts 
with a big enough graph and removes as many elements as possible. However, for edges, 
besides the actions of the productions on edges we need to keep track of the actions of 
the productions on nodes because some potential dangling edges may become available 
(if their incident nodes are added by some grammar rule). 

Notice that Th. 4.4.2 in [15] proves that the certainty part of the initial digraph is 
the one that appears in eq. ( [33] ). For the nihil term X/iCxTxKy it is easier to think in 
what must be or must not be found in G. 

We proceed by induction on the number of productions. For the time being, for 
simplicity, we omit the effect of adding nodes which may turn potential dangling edges 
into available ones. In a sequence with a single production it should be obvious that Ki 
(and only Ki) needs to be demanded. 

For a sequence of two productions §2 = P2',Pi, Ki is again necessary. It is clear 
that Ki V K2 with KY KY = - i.e. all nodes and hence edges unrelated - would 
be enough, but it may include more elements than strictly needed. Among them, those 
already deleted by pi (once they are deleted they belong to G) and those that already 
appear in Ki and that are not added by pi - fT A'l -. If these elements of K2 are not 
going to be considered, we need to and their negation: ei {riKi)K2- Altogether we get 
Ki V ei(jiKi)K2- Some simple manipulations prove that: 

Ki V K2ei{riKi) = Ki v K2ei (ri v Ki) = 

= Kiv K2(ei wTiKi) = KYQ^ . (34) 

Minimality is inferred by construction. If any other element was removed then either 
Pi^ or could not be applied (and still consider dangling edges). It is not difficult to 
check that the sequence P2^',Pi'^ can be applied to Ki v K2Qi. The expressions for 



sequences of three, four, n poductions are: 



N3 = N2V XgraQi Q2 



N4 = N3V K4r3r2Qi raQa Q3 

n 

n 



i=2 



(35) 
(36) 

(37) 
(38) 



There are two tricky steps. The first one is how to derive Nn in eq. dSTT i and the 
second is how to obtain its equivalent expression dSST l. The reader is referred again to 
the aforementioned proof in Iil5il where detailed explanations are given for a similar 
case. 

Once we get here it is easy to obtain Vi ^xKy). First, note that the sequence is 
coherent so the third term in eq. ( |38] | is zero. Second, as Ki = i^i v ri, the ri can be 
simplified because a v a6 = a v 6 in propositional logic. 

Finally, the same reasoning applies for those nodes that are added. So we do not 
only need to remove elements erased by previous productions but also edges that are 
not incident to any non-existent edge, v" {^xKy) ^ Vi {^xTxKy)M 
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Fig. 11. sequence of Two Productions 



Example. oFigurelTTI includes two productions with their nihilation matrices Ki and 
K2 ■ The initial digraph of the sequence s = 52 ; 91 is 



M{s2) = Vi (r^LyviexT^Ky) = (riLi vrir2£2) vi (eiTiKiv eie2TiT2K2) 



{Li V riL2)vi (TiKi v eiTiT2K2 
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We have represented Mc{s2) v Mm{s2) to the left of Fig. fT2] together with its 
evolution as well as the final state, S2 (Af (52)). To the right of the same figure there is 
the same evolution but limited to edges and from the point of view of swaps. With black 
solid line we have represented the edges that are present and with red dotted line those 
that are absent. Recall that swaps interchange them.B 
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Fig. 12. Initial Digraph of a Sequence of Two Productions Together with its Evolution 



A final remark is that T makes the number of edges in G as small as possible. For 
example, in eie2TiT2K2 we are in particular demanding eiTiT2r2 (because K2 = 
r2 V €202). If we start with a compatible host graph, it is not necessary to ask for 
the absence of edges incident to nodes that are added by a production (we called them 
potentially available above). Notice that these edges could not be in the host graph as 
they would be dangling edges or we would be adding an already existent node). 



7 Compatibility and Congruence 

This section revises some more sequential results, adapting and extending them via 
MCL. The notions we cope with are the image of a sequence, compatibilit}0 and G- 

^' Compatibility for a single production has been tackled in Prop.|4l Sec. ID 



congruence. By the end of the section we will very briefly touch on sequential indepen- 
dence, application conditions and graph constraints. 

The image of a sequence of productions s = p„; . . . ; pi acting on its initial digraph 

M(s) = Mc(s) V iMN(s) is given by: 

s (M(s)) = (^A'^e^ry v /\ejMcis)^ v z (^A^^r^Cy v /\ rjMAr(s)^ . (39) 

Equation ( [39] ) is deduced by simply applying each production to the initial digraph. 
We would like to interpret the shape of the image of a sequence as a production by 
setting s = e(s) v ir(s), with e(s) = Vr=i ^^'^ ''^i^) ~ '^i e^^'yEl We could then 
put it as (A/(s), s) and for example calculate its associated swap: 

n n 

w{s) = /\ (ejfj) V i Y (ej v r^) . 

Unfortunately this is not possible because although in the certainty part of eq. ( [39] l 

r(s) = AiCxTy and e(s) = \/r=i '^^e nihil part we find Aj=i ^ ''(■s) 
AiTxCy e(s). 

Compatibility asks for "closedness" of the space (graphs) with respect to the speci- 
fied operations. In essence, it demands the lack of dangling edges. Definitions of com- 
patibility for increasingly general concepts can be found in fTSl: single simple digraph, 
production and sequence. According to Prop. |4] productions act on edges and on ver- 
tices. They are obviously related but this relation has not been demonstrated. It is of 
importance in order to study the evolution of the nihil part of complex terms. What one 
production forbids, another production may need or even can make accessible again. 

Proposition 5 (Compatibility). Let s = pn', ■ ■ ■ ] pi be a sequence made up of compat- 
ible productions. If 

V? (e.r.Mc(s,)Mw(s,)) = (40) 

then s is compatible, where Mc{sm) and MN{sm) are the certainty and nihil parts of 
the initial digraphs of Sm = Pm', ■ ■ ■ ',Pi, me {!,..., n}. 

Proof (Sketch) 

□Equation ( l40l l is a restatement of the definition of compatibility for a sequence of 
productions. The condition appears when the certainty and nihil parts are demanded to 
have no common elements. Compatibility of each production is used to simplify terms 
of the form LiKiM 

Compatibility and coherence are related notions but only to some extent. Coher- 
ence deals with actions of productions while compatibility with potential presence or 
abscense of elements. This is better understood if we think in derivations: when the left 
hand side L v iK of rule p is matched in a host graph G v iG, all elements of L must 
be found in G and all edges of K must be found in G. When p is applied a new graph 



The idea behind would be the composition of a sequence of productions to derive a single 
production: s„ = pn; ■ • • ; pi i-^ c„ = p„ o . . . o pi. 



H V iH is derived. Again, all elements of R have to be found in H and all edges in Q 
will be in H, no matter if some of them are now potentially usable (say p adds some 
nodes and some potentially dangling edges are not dangling edges anymore). 

Now we turn to G-congruence, which studies equality of initial digraphs for a 
sequence s = p„; . . . ;pi and a permutation of it, s' = (j{s). All the job for ad- 
vancement and delaying of productions - permutations and 5, wherq3 advance- 
ment is (/) = (1 2 ... n — 1 n) and delaying is (5 = (n n — 1 ... 2 1), i.e 
(/)(s) = p„_i;p„_2 . . . ;pi;p„ and 8{s) = • . .P2 - is done in [151, Sec. 6.1, 

so we state the result without proof. 

Theorem 4 (G-congruence). With notation as above, sequences s and (j){s) are G- 

congruents if F = v iF~ = 0, where 

F+ = LnVl-^e^Ky {ry v e„) and F- = RnV'^-'^r^Ly {cy v r„) . (41) 
Also, s and 6{s) are G-congruents if D = v iD~ = 0, with 

D+ = LiS/l^e^Ky (ry v ei) and = KiV^r^Ly {cy v ri) . (42) 

Proof 
□■ 

An easy remark is that the complex term C+ v iC^ in Th.|2]provides more infor- 
mation than just settling coherence as it measures non-coherence: Problematic elements 
(i.e. those that prevent coherence) would appear as ones and the rest as zeros. The same 
holds for F+ viF~ and Z?+ viD~ in Th.|4]for congruence and eq. (l40b in Prop.|5]for 
compatibility. 

There are some relevant topics that we have not mentioned such as sequential in- 
dependence, application conditions and graph constraints. We briefly discuss how they 
could be handled wifli MCL. 

With respect to the image of a sequence and sequential independence, recall that 
MCL naturally uses swaps rather than productions. This abstraction has its effects on the 
interpretation of operations. On the positive side, among many other things, swaps are 
a nice redefinition and generalization of productions that take into account the certainty 
and nihil parts; on the negative side, our intuition needs to be adjusted. For example, 
consider a production p that only deletes edge (1, 2) and does nothing else. Suppose 
that it is applied twice to the graph G that consists of nodes 1, 2 and edge (1, 2). In this 
case p;p(G) = G which is algebraically correct. However, it does not encode "delete 
edge (1, 2) twice". Of course, the point here is that of completion: we would rather 
have considered its application to G', made up of nodes 1, 1' and 2 and edges (1, 2) and 
(1', 2). A similar reasoning shows that sequential independence is "granted" if we rely 
only on algebraic operations and do not pay attention to completion: 

P2;Pi{C) = ((C,P{pi)),P{p2)) = CP{pi)P{p2) = CP{p2)Pipi)=pi;p2{C). 



Numbers in the permutation refers to the position that the production occupies inside the se- 
quence, not to its subindex. 



Previous comments highlight some of the reasons why coherence, compatibiHty, initial 
digraph and G-congruence are so valuable, justifying their inclusion and also linking 
the present and previous sections to Sec.|5] 

Regarding application conditions and graph constraints, they are not difficulty re- 
lated to what has been presented so far. Recall from Sec. |5] that swaps transform ele- 
ments in the same diagonal of the Sierpinski gasket. If they are allowed to be applied to 
g € <3 instead of the restricted case that we have studied (Sj), we may impose limits on 
what elements can not be added nor deleted by sequences of productions (swaps). This 
is because if one edge is in the certainty part and in the nihil part, it can not be deleted 
by any swap. On the contrary, if one edge does not appear neither in the certainty nor 
in the nihil parts, it is not possible for a swap to add it. 

If we call any of these situations a swap restriction, it can be guaranteed that a 
sequence will not add nor delete (or both) some element, despite the actual definition 
of the productions that make up the sequence or the grammar Again, geometrically, we 
are choosing the diagonal inside the Sierpinski gasket in which all operations will take 
place. 

8 Conclusions and Future Work 

In this paper we have introduced Monotone Complex Logic (MCL) which, in our opin- 
ion, is an interesting topic in itself. With respect to Matrix Graph Grammars (MGGs), 
MCL allows the encoding of simple digraphs and grammar rules using complex terms. 
We believe it is a natural representation in the MGGs context, as productions act on pairs 
of graphs (L, K) ^ (p(L),p^^(K)). Relevant algebraic structures for their study have 
been introduced (PMCA, PMMA, ®, Sj). Swaps allow studying and classifying pro- 
ductions according to their dynamic behaviour, defining a surjective morphism into the 
self-adjoint graphs in Sj. 

The rational enconding of Boolean complexes gives an embedding of MGGs into 
a subset of the complex numbers: the Sierpinski gasket. Using such representation we 
have been able to introduce standard geometric and analytic concepts such as a scalar 
product, a norm and a notion of distance in MGGs. This generalizes the theory de- 
veloped in fTSl and opens the door to the study of dynamics of infinite graphs with a 
countable number of nodes (this topic is left for future research). Finally, some of the 
most relevant concepts of MGGs have been expressed and reinterpreted using MCL: 
coherence, initial digraph, image of a sequence, compatibility and G-congruence. 

Our main interest is complexity theory so we have to introduce a measure of com- 
plexity. The natural proposal seems to be the geodesic distance, which measures the 
cost of reaching one element from another one through elements of the MGG (rules of 
the grammar). However, the natural distance here is not the Euclidean one, but the xor 
metric restricted by available operations. 

One foreseen advantage of the results in this paper is that there is a lot of interest 
and current research activity on the Sierpinski gasket 12171181 . In the mid-long term we 
plan to continue our work towards computational complexity theory through MGGs. 
It is our opinion that one of the main "problems" of current approaches to complexity 
theory is that there are very few links to other branches of mathematics (there are some 



exceptions though, such as IfTOl ). As MGGs are a compact and path connected fractal, 
it seems promising to introduce harmonic and functional analysis and noncommutative 
geometry techniques, apart from those already available in MGGs. Two main research 
directions will be explored in the future: measurable Riemannian geometry as in \J] 
and noncommutative geometry as in ||2l • 

Notice that it is not difficult to interpret MGGs as a model of computation (we 
are preparing a paper on this topic). Also, it might be of interest to encode properties 
of graphs (such as coloring) using graph grammars, translating static properties into 
equivalent dynamic properties of associated sequences. 

Another point of interest might be the introduction of stochastic analysis. This is 
closely related to MGGs as a model of computation and the way grammar rules are 
selected (a source of non-determinism). Other source of non-determinism appears in 
case there are several places in a host graph to which a production can be applied. 

There are many more topics for further research, e.g. graph constraints, derivations, 
applicability, reachability, dynamic encoding of static properties and infinite graphs 
some of which we have already commented on. 
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